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1.  Introduction 

Das  (1956)  presents  a  method  of  evaluating  the  integral 
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(where  if ( x^  ,  Xg  .  •••  ,  xn)  is  the  joint  multivariate  normal  density  function 
with  zero  means  and  nonsingular  variance-covariance  matrix  £)  through  the 
combining  of  n  +  k  independent  normal  variables  with  zero  means  and  unit 
variances.  Later  Marsaglia  (1963)  shows  that  this  is  a  special  case  of  a 
convolution  formula.  The  complexity  of  implementing  the  solution  is  highly 
dependent  upon  the  size  of  k  and  Marsaglia  (1963)  notes  that  k  equal  to 
n  minus  the  multiplicity  of  the  smallest  latent  root  of  £,  can  always  be 
achieved.  This  note  investigates  properties  of  fc  that  will  allow  smaller 
values  of  k  . 


2.  The  equivalent  expression 

As  a  method  for  evaluating  I  ,  Das  (1956)  considers  two  row  rectors 
y'  =  (Yj  .  •  yn>  and  z'  =  (Zj  •  z2  •  ***  •  whose  elements 

are  normally  and  independently  distributed  with  zero  means  and  unit  variances. 
The  problem  is  to  choose  a  positive  constant  c  and  an  n  x  k  real  matrix  B 
such  that 
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r  n 


i) 


2 


for  then  x*  =  (Xj  ,  xg  ,  ••  •  .  xR)  can  be  expressed  as 

x  =  cy  -  Bz  . 

I  can  now  be  expressed  as 
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where  P[a]  =  Pr(yj  s  a)  .  It  is  evident  that  a  small  k  is  advantageous. 
The  slight  change  in  1) 


t  =  +  BB’  , 


3) 


where  C  is  a  diagonal  matrix  with  positive  diagonal  elements  c.  ,  results 
in 

x  =  Cy  -  B z 

and  2)  of  the  form 
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3.  Latent  vector  properties 

If  £  can  be  expressed  in  the  form  3)  and  w^  is  a  latent  vector 
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of  £  then 
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where  bj  is  the  jth  column  of  B  .  From  expression  5) 
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where  w.  is  the  ith  element  of  w  and  0.  =  b|w  .  If  the  c?  are 

not  all  equal  the  n  equations  of  6)  are  parametric  equations  for  a 

p 

k  dimension  variety  in  n  space  of  order  two,  ;  a  brief  discussion 

of  this  may  be  found  in  Kendall  (1956),  page  5.  That  is  given  a  vector 

w^  ,  algebraic  operations  to  eliminate  and  the  0^  will  result  in 

o 

n-k  second  degree  homogeneous  equations  in  the  w..  .  If  the  c.  are 
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all  equal  say  to  c  then  6)  can  be  written  as 
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where  =  kjw^Ax^  -  c  )  ,  v:hi/*h  are  the  parametric  equations  of  a 
k  dimension  flat,  . 

Conversely  if  the  latent  vectors  of  £  all  fall  in  a  of  the  form 

of  6)  with 


4  0  *nd  ci >0 


then  £  can  be  written  in  the  form  3)  .  If  the  latent  vectors  of  £  all 
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fall  in  a  with 


1 0  ’> 

then  it  can  be  written  in  the  form  1)  . 

Thus  the  smallest  value  of  k  for  evaluating  I  through  form  4)  is 

2 

determined  by  the  smaller  of  the  minimum  dimension  of  the  Vk  ’s  of  the 
form  6)  satisfying  conditions  8)  and  the  minimum  dimension  of  the 
satisfying  the  condition  9)  that  contain  the  latent  vectors  of  £  . 

Alternative  conditions  can  be  given  for  the  special  case  of  k  =  1  . 

6)  becomes 


or 
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10) 


which  is  the  parametric  expression  of  a  plane,  S2  ,  passing  through  the 
origin.  Thus  if  the  s'^  =  (l/w^  ,  i/*^  •  *’*  •  l/*»  all  fall  in  a 
plane  with  parametric  conditions  10)  satisfying 
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2  2 

only  one  auxiliary  variable  is  necessary.  For  the  case  of  c.  =  c  for  all  i  , 

2 

n-1  of  the  equal  c  and  their  corresponding  0^  are  zero.  That 
is  the  minimum  latent  root  of  £  has  multiplicity  k-1 

4.  Summary 

In  general  the  minimum  number  of  auxiliary  variables  necessary  to 


evaluate  a  multivariate  normal  integral  through  the  method  of  Das  (1956) 


5 


can  be  determined  by  the  dimensions  of  the  second  order  varieties  with  para¬ 
metric  form  6)  and  the  dimensions  of  the  flats  that  contain  the  latent 
vectors  of  the  variance-covariance  matrix. 

This  problem  could  also  be  looked  at  as  determining  the  diagonal 
matrix  D  such  that  d£d  has  maximum  multiplicity  of  its  smallest 

latent  root.  Then  Marsaglia' s  (1963)  solution  suffices.  It  is  interesting 
to  note  that  the  difficulty  of  the  solution,  the  dimension  of  the  integration 
in  I  ,  is  affected  by  change  of  scales  of  the  normal  variables. 
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zero  means  and  nonsingular  variance-covariance  matrix  $)  through  the  combining  of 
n  +  k  independent  normal  variables  with  zero  means  and  unit  variances.  Later 
Marsaglia  (1963)  shows  that  this  is  a  special  case  of  a  convolution  formula.  The 
complexity  of  implementing  the  solution  is  highly  dependent  upon  the  size  of  k  and 
Marsaglia  (1963)  notes  that  k  equal  to  n  minus  the  multiplicity  of  the  smallest  latent 
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allow  smaller  values  of  k. 
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